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In this article, we analyze some nuclear structure aspects of the Ou double beta decay nuclear 
matrix elements (NME). We give results for the decays of ''**Ca, '^'^Ge, **^Se, ^^^Sn, ^^*Te, "°Te, and 
^•^^Xe, using improved effective interactions and valence spaces. We examine the dependence of the 
NME's on the effective interaction and the valence space, and analyze the effects of the short range 
correlations and the finite size of the nucleon. Finally we study the influence of the deformation on 
the values of the NME's. 
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I. INTRODUCTION 

The double beta decay is the rarest nuclear weak pro- 
cess. It takes place between two even-even isobars, when 
the decay to the intermediate nucleus is energetically for- 
bidden or hindered by the large spin difference between 
the parent ground state and the available states in the 
intermediate nuclei. It comes in three forms: The two- 
neutrino decay /3/92i/ is just a second order process me- 
diated by the weak interaction. It conserves the lepton 
number and has been already observed in a few nuclei. 
The second mode, the neutrinoless decay f3l3o,j, needs an 
extension of the standard model of the electroweak in- 
teractions as it violates lepton number. A third mode, 
/3/3oiy,x is also possible in some extensions of the standard 
model and proceeds via emission of a light neutral boson, 
a Majoron x- The last two modes, not yet experimentally 
observed, require massive neutrinos -an issue already set- 
tled by the recent measures by Super-Kamiokande ^ll] , 
SNO Hi and KamLAND Interestingly, the double 
beta decay without emission of neutrinos would be the 
only way to sign the Majorana character of the neutrino 
and to distinguish between the different scenarios for the 
neutrino mass differences. In what follows we shall con- 
centrate mostly in the PPou mode. The uprising of inter- 
est in the observation of the PPoi, decay was somewhat 
obscured by the analysis of the status of the calculations 
of the nuclear matrix elements, that enter in the lifetime 
of the decay together with the effective neutrino mass, 
made in ref. [H. The authors concluded on a very pes- 
simistic note, saying that the spread of the available val- 
ues for the nuclear matrix elements was such that there 
was no hope to translate the experimental signal into use- 
ful input for the theories beyond the standard model. A 
critical assessment of the results of the many calculations 
available in the framework of the quasi particle random 
phase approximation (QRPA) has been made recently 
[5|, with a much more optimistic conclusion. In this ar- 
ticle we want to continue improving upon the reliability 
of the calculated nuclear matrix elements (NME) in the 
framework of large scale applications of the Interacting 
Shell Model (ISM) . In addition, we put forward a process 
of benchmarking the different approaches, and study the 



stability of the NME's under reasonable modifications of 
the nuclear structure inputs of the calculations. 



II. DOUBLE BETA DECAYS 

As we have already mentioned, some nuclei, otherwise 
nearly stable, decay emitting two electrons and two neu- 
trinos {2v f3f}) by a second order process mediated by 
the weak interaction, that has been experimentally mea- 
sured in several favorable cases. The decay probability 
contains a phase space factor and the square of a nuclear 
matrix element 



(1) 



If the neutrinos are massive Majorana particles, the 
double beta decay can take place without emission of 
neutrinos {Ov P(3) . In this case the transition is mediated 
by terms that go beyond the standard model. The decay 
probability contains a phase space factor ^ , the effective 
electron neutrino mass (a linear combination of the mass 
eigenvalues whose coefficients are elements of the mixing 
matrix) and the nuclear matrix element 0, Q • 
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The 21/ matrix element can be written as: 



{GDiJ)\at-\lt){lt\at-\F) 
AE,, 



(3) 



The nuclear structure information entering in this ma- 
trix element consists on the wave function of the ground 
state of the father nucleus (F) and the wave function of 
the state of the grand daughter, GD(J), nucleus (J=0, 
J=l, or 3—2) to which the decay proceeds. In addi- 
tion, the wave functions and excitation energies of all 
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the l"*" states in the odd-odd daughter nucleus are in 
principle necessary. The spin-isospin operators in the 
nuclear medium are quenched by a factor q«l/g^, and 
this quenching factor is also required to reproduce the 
experimental data of the 2v double beta decays. 

To compute exactly the {Sv matrix element, we would 
also need the wave function of the ground state of the 
father nucleus and of the ground state grand daughter, 
but in addition we would need all the wave functions 
and excitation energies of all the J'^ states in the odd- 
odd daughter nucleus. Most conveniently, the matrix 
elements can be approximately obtained in the closure 
approximation, that is good to better than 90% due to 
the high momentum of the virtual neutrino in the nucleus 
(wlOO MeV) Hence, the matrix element can be written 
as: 



82, 



Se 



0* 1* 2* 3* A* 5* 6* 1* 
1 2 3 4 5 6 7 



= {GD\h(\r\ - r2|)(ai • o^m^)\F) (4) 

where r^j) is the neutrino potential («l/r). In this 

approximation no knowledge of the intermediate nucleus 
is directly implied. 

The transition operators are usually obtained from the 
Hamiltonian of Doi et al. [il However, according to recent 
claims [oj, additional terms originating in the coupling 
to the virtual pions, should give non-negligible contribu- 
tions. The finite size of the nucleon and the short range 
correlations need to be taken into account in the calcula- 
tion of the two body matrix elements of the Oi^ two-body 
transition operators as well. 

Perhaps the most relevant issue for our mastering on 
the neutrinoless double beta decays is to get a better 
insight in the physical content of the two body transition 
operator. Only with this knowledge could one decide 
whether the nuclear wave functions that enter into the 
calculation of the NME's contain the degrees of freedom 
that correspond to the correlations to which the operator 
couples dominantly. 

The two body decay operators can be written generi- 
cally as: 

Mi^^) (5) 

In words, what the operator does is to annihilate two 
neutrons in the parent nucleus and to create two protons 
The NME is the overlap of the resulting object with the 
grand daughter ground state. 

The contributions to the Ov matrix element as a func- 
tion of the J of the of the decaying pair have a very telling 
structure as can be seen in Figure [TJ The dominant con- 
tribution corresponds to the J=0^ pairs, while all the 
other pairs have much smaller contributions, but all of 
them of sign opposite to the leading term. If the ini- 
tial and final wave functions had seniority zero, only the 
leading term would contribute and the matrix element 
would be maximal. This is a first indication, relating 



FIG. 1: The matrix element of the * Se Kr decay as a 
function of the J of the decaying pair 



the pairing content of the nuclear wave functions and the 
neutrinoless double beta decay operator. 



III. INTERACTING SHELL MODEL (ISM) AND 
QRPA CALCULATIONS 

In the quest for better wave functions to describe the 
double beta decay processes, two main avenues have been 
explored. The interacting shell model in larger and larger 
valence spaces, with ever improving effective interactions, 
and the quasiparticle RPA. Ideally, both methods should 
be able to produce good spectroscopy for parent, daugh- 
ter and grand-daughter, even better if it extends to a full 
mass region, correct total Gamow Teller strengths and 
strength functions, 2^ matrix elements, etc. In brief, the 
goal is to have a description as close to perfect as possible 
of the dynamics of the nuclei involved in the transition. 
Until rather recently, ISM calculations that could encom- 
pass most of the relevant degrees of freedom of the nuclei 
of interest, were only available for the lighter emitters, 
and, as a consequence, most systematic studies were per- 
formed in the QRPA framework. As we will discuss in 
the next section, all the potential double beta emitters 
are now within the reach of the ISM except ^^"Nd, that, 
being deformed, is also out of the reach of the spherical 
QRPA calculations. 

In the ISM approach, the valence spaces contain a 
number of orbits that is "small" compared to the QRPA, 
however, all the possible ways of distributing the valence 
particles among the valence orbits are taken into account. 
In the QRPA calculations, the number of active orbits is 
larger, but only Ip-lh and 2p-2h excitations from the 
normal filling are considered (and not all of them) . 

The effective interactions used in the ISM calculations 
are usually G-matrices whose monopole behavior is fitted 
to the spectroscopic properties of a large region of nuclei. 
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in general those comprised between to magic closures for 
the neutrons and for the protons. In some cases the inter- 
actions are plainly fitted to a set of experimental masses 
and excitation energies. In the QRPA description, the 
starting point is provided by realistic or schematic in- 
teractions, however, both the particle-particle and the 
particle-hole channels of the interaction are affected of 
strength parameters dubbed gph and gpp that are fitted 
to selected experimental data. In particular, it has been 
shown that the double beta decay matrix elements de- 
pend critically of gpp . There are different opinions among 
the QRPA practitioners about the best choice of the value 
of gpp. The two preferred options being to fit the exper- 
imental 2i/ matrix elements, or to fit the Gamow Teller 
strength functions. In a sense, none is really safe, be- 
cause the pure at^ channel plays a rather minor role in 
the Ov decay. 

The dominant correlations in the nuclei are due to 
the pairing and quadrupole-quadrupole terms of the nu- 
cleon nucleon interaction in the medium. In the ISM, 
the pairing correlations are treated exactly within the 
valence space. Proton and neutron numbers are exactly 
conserved. Proton-proton, neutron-neutron, and proton- 
neutron (iso-vector and iso-scalar) pairing channels are 
included on equal footing. On the contrary, in the spher- 
ical QRPA only proton-proton and neutron-neutron pair- 
ing terms are considered. They are treated in the BCS 
approximation. Proton and neutron numbers are not ex- 
actly conserved. From the spherical shell model point 
of view, this approach is a seniority zero approximation 
at the BCS level. When the RPA correlations are taken 
into account, higher seniority components appear in the 
ground states of the father and grand daughter nuclei, 
however, for the RPA to be valid, their amplitudes must 
decrease rapidly with increasing seniority. 

The multipole correlations and the eventual intrinsic 
deformation can be properly described in the laboratory 
frame in the ISM description. Angular momentum con- 
servation is exactly preserved. In the QRPA the multi- 
pole correlations of the ground state are treated at the 
RPA level. The posibility of having permanent intrinsic 
deformation is not contemplated yet. 



IV. THE ISM VALENCE SPACES RELEVANT 
FOR fjp DECAYS 

Classical Ohuj ISM valence spaces are the p, sd and 
pf shells. These correspond to the harmonic oscillator 
major shells of principal quantum number p equal to 
1, 2, and 3. We use the following notation; in shell p, 
the orbit j=p+l/2 is denoted by the spectroscopic label 
of l=p and the remaining ones by rp. The space rag is 
adequate for ^^Ge, ®^Se and their descendants. With the 
r4h space for the neutrons and the g-r4 for the protons, 
we aim to understand the decays of ^^Zr, ^°"Mo,^-'^°Pd, 
and ^-'^^Cd. Finally the r4-h space is a natural choice for 
the description of the ^^"^Sn, ^^^~^^°Te, and ^■^^Xe decays. 





TABLE I: 


: Update of the ISM Ou results 




(nv) 


for T 1 = 

2 


10"" y. 






128 rpg 

130,pg 

i3«Xe 


0.85 
0.90 
0.42 
0.45 
1.92 
0.35 
0.41 




0.67 
2.35 
2.26 
2.11 
2.36 
2.13 
1.77 


1.14 
1.10 
1.10 
1.13 
1.13 
1.13 
1.13 



All these spaces are accessible to larg e scale ISM 
descriptions. The Strasbourg-Madrid [1^ Shell-Model 
codes can deal with problems involving bases of O(10^°) 
Slater determinants, using relatively modest computa- 
tional resources. 



V. UPDATE OF THE ISM Oi^ RESULTS 

In the valence spaces Ts-gg/2 (^^Ge, ^^Se) and r4-/iii/2 
(i24Sn, i28-i30rpg^ i36xe) we have obtained high quahty 
effective interactions by carrying out multi-parametrical 
fits [ll| whose starting point is given by realistic G- 
matrices [l^]. With these interactions the beta decay 
properties of a large set of nuclei are well reproduced. 
The 2i/ double beta decay half-lives are found in rea- 
sonably good agreement with the experimental results 
as well. In the valence space proposed for ^^Zr, ^"^Mo, 
^^°Pd and ^^^Cd, our results are still preliminary and 
subject to further improvement both on the interaction 
side and on the removal of the yet necessary truncations. 

Our results are obtained in the closure approximation, 
with the short range and finite size corrections mod- 
eled as described in [l^; using ro=1.2 fm to make the 
matrix element dimensionless; with gyi=1.25; and with- 
out higher order contributions to the nuclear current. A 
preliminary estimation of the higher order contributions 
gives a reduction of the ISM NME's in the range of 10% 
. Our present best values are collected in table [B xf is 
defined as: 

X.= (g)'Mr (6) 

Except in the case of doubly magic **Ca, whose NME 
is severely quenched, all the other values cluster around 
a value M'^°'')=2.5. The limits of the effective neutrino 
mass for a half life limit of 10'^^ y, that incorporate the 
phase space factors, show a mild preference for some of 
the potential emitters. 

In table |TT] we compare the ISM results with the most 
recent QRPA calculations including the higher order cor- 
rections discussed before. The range of values of the 
NME's shown in the table is that given by the authors. 
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TABLE II: Comparison of the present ISM results with the 
range of values proposed in different QRPA calculations; (a) 
ref. [^,[1^, (b) ref [l^. The values of column (a) have been in- 
creased by 10% with respect to the numbers given in ref. [T4| . 
for a proper comparison with the others, because of their use 
of ro=l.l fm. In the same vein, the ISM values should also be 
reduced by ^10% to account approximately for the absence 
of higher order corrections to the nuclear current. 



Mo. ISM QRPA(a) QRPA(b) 



^«Ge 


2.58 


3.81 - 


4.96 


4.06 - 


5.26 


»^Se 


2.49 


3.20 - 


4.42 


2.72 - 


3.60 


128 


2.67 


2.79 - 


4.00 


3.37 - 


4.42 


130,pg 


2.41 


2.57 - 


3.59 


3.00 - 


4.22 


i3«Xe 


2.00 


1.39 - 


2.32 


1.92 - 


2.92 



and derives from the different choices of gpp and gA used 
in the calculations, as weU as from the use or not of a 
renormahzed version of the QRPA. In addition, in order 
to make the comparison more transparent, we have se- 
lected the results that treat the short range correlations 
by means of a Jastrow factor. Overall, the two sets of 
QRPA calculations are now compatible. The ISM pre- 
dictions of the NME's are systematically smaller than 
the QRPA central values, except in the case of ^"^^Xe. 
This nucleus is semi-magic, and the present experimen- 
tal limit on its 21/ decay half life is surprisingly large, 
perhaps indicating that some subtle cancelation mecha- 
nism is at work. For the others, a plausible explanation 
of the discrepancy, relating it to the implicit seniority 
truncations present in the spherical QRPA calculations, 
will be presented elsewhere. 



VI. EXPLORING THE DEPENDENCES OF 
THE ISM RESULTS 

A. Effective interaction 

The ISM results depend only weakly on the effective 
interactions provided they are compatible with the spec- 
troscopy of the region. For instance, in the pf shell we 
have three interactions that work properly, KB3 p^ . 
FPD6 [13 and GXPFl Their predictions for the 

21/ and the neutrinoless modes are quite close to each 
other (see table IIII| ). Similarly, in the r^g and r4h spaces, 
the variations among the predictions of spectroscopically 
tested interactions are small (10-20%). 



B. Finite size and short range corrections 

There has been a certain debate among QRPA practi- 
tioners about the amount of the reduction of the NME's 
due to the short range correlations ,19]. The debate has 
become milder after the publication of ref. [l^. In the 



TABLE III; Comparison of the ISM results for the ***Ca de- 
cays using the effective interactions KB3, FPD6, and GXPFl 





KB3 


FPD6 


GXPFl 


Mgt{2u) 


0.08 


0.10 


0.11 


Ugt{Ou) 


0.67 


0.73 


0.62 



TABLE IV: Influence of the short range correlations (SR) and 
of the nucleon finite size (FS) on the nuclear matrix elements 
of the neutrinoless double beta decays 



Ugt{Ou) bare +FS +SR +FS+SR 



^«Ca ^*«Ti 


1.21 


0.93 


0.74 


0.67 




3.55 


2.88 


2.57 


2.35 


82 Se ^«2Kr 


3.33 


2.72 


2.43 


2.26 


i30xe ^i30Xe 


3.34 


2.70 


2.32 


2.13 


i36xe ^i36Ba 


2.76 


2.24 


1.93 


1.77 


ISM description, 


once 


the finite 


size of the 


nucleon is 



taken into account by means of the standard dipole form 
factor |20|], the effect of the short range correlations, that 
we model by a Jastrow ansatz following ref. [2l[ , is about 
one half of what it would be without it. In fact, the 
ISM corrections for finite size and short range effects 
with the Jastrow prescription are quite close in per- 
centage to the QRPA values. Short range and finite size 
corrections proceed mainly through the reduction of the 
J'^^O^ pair contribution. Preliminary ISM calculations 
using a softer prescription for the short range correlations 
do not change this picture at all. In fact, the relative val- 
ues of the ISM NME's with respect to the QRPA ones 
seem to be independent of the choice of the (common) 
prescription used for the the short range correlations. 

C. The influence of deformation 

Changing adequately the effective interaction used in 
the calculations, we can increase or decrease the deforma- 
tion of parent, grand-daughter, or both, at will. In this 
manner, we can gauge the effect of these variations on the 
decays. We have artificially changed the deformation of 
^^Ti and "'^Cr adding an extra XQ Q term to the effective 
interaction. The results are presented in Fig. [51 Positive 
values of A increase the deformation while negative values 
reduce it. For zero values of both A's, ^*Cr is already well 
deformed, while "^^Ti is transitional. The circle to the up- 
per left corresponds to the spherical-spherical situation, 
the square to the upper right to equally deformed Tita- 
nium and Chromium, and the diamond to the bottom left 
to a spherical Titanium and a very deformed Chromium. 
Our conclusion is that a large mismatch of deformation 
can reduce the /3/3 matrix elements by factors as large 
as 2-3. This exercise indicates that the effect of defor- 



5 




-2-10 1 



FIG. 2: Influence of deformation in the (flctitious) decay of 

mation is very important and cannot be overlooked. We 
have reached similar conclusions for heavier emitters and 
a systematic exploration in realistic cases is under way. 



stant, even if the total Gamow- Teller strengths, (GT+) 
and (GT-), increase from 0.15 to 0.34 and from 20.5 to 
26.9. The larger GT+ strength is somehow compensated 
by a larger cancellation among the contributions of the 
different intermediate states and by a shift of the centroi'd 
of the strength toward higher energy. 

We have also computed the ^^^Xe decay in the r4h 
space including 2p-2h excitations from the Og| proton or- 
bit and the matrix element increases less than 10%. In 
another set of calculations, we have included 2p2h neu- 
tron excitations toward the Oh 9 and If?. The occupan- 
cies that we obtain are relatively large (0.25 neutrons in 
each orbit) and the effect is to increase the matrix ele- 
ment by 15%. It is interesting to note that the increase 
with the two orbits simultaneously active is equivalent to 
that obtained including one or another orbit separately. 
Therefore there is no pile-up of the contributions of the 
small components of the wave function. As a preliminary 
conclusion, the ISM results seem to be robust against the 
inclusion of small components of the wave function. 



VIII. CONCLUSIONS 



VII. ISM CALCULATIONS IN QRPA-LIKE 
VALENCE SPACES 

The ISM valence space for the ^^Ge and ®^Se decays 
has traditionally encompassed the orbits Ipa, Ofs, Ipi, 
and Og£. In the QRPA, two major oscillator shells are 
taken into account; Of?, Ipa, Ofs, Ipi, Ogg, Ids, Og?, 
2s 1, and Ida. 

2 2 

As a first step toward a more complete benchmark- 
ing, we have evaluated the influence of the 2p-2h jumps 
from the lf7/2 orbit - ^^Ni core excitations- in our results 
for the *^Se decay. Similar calculations for the ^^Ge de- 
cay are under way. The calculation in the full rag space 
plus 2p-2h proton excitations from the Of 7 orbit gives a 
20% increase of M^'^ , but probably we overestimate the 
amount of core excitations. Our Ofi proton occupancies, 
7.71 and 7.69 in ^^Se and ^^Kr are smaller than the BCS 
occupancies of Rodin et al. Q, 7.84 and 7.84. There- 
fore the above 20% must be taken as a very conservative 
upper bound. The 21/ matrix element remains nearly con- 



Large scale shell model calculations with high quality 
effective interactions are available or will be in the imme- 
diate future for all but one of the neutrinoless double beta 
emitters. The theoretical spread of the values of the nu- 
clear matrix elements entering in the lifetime calculations 
is greatly reduced if the ingredients of each calculation 
are examined critically and only those fulfilling a set of 
quality criteria are retained. A concerted effort of bench- 
marking between ISM and QRPA practitioners would be 
of utmost importance to increase the reliability and pre- 
cision of the nuclear structure input for the double beta 
decay processes. 
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